Kondo spin liquid and magnetically long-range ordered states in the Kondo necklace 

model 
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A simplified version of tlie symmetric Kondo lattice model, 
the Kondo necklace model, is studied by using a representa- 
tion of impurity and conduction electron spins in terms of lo- 
cal Kondo singlet and triplet operators. Within a mean field 
theory, a spin gap always appears in the spin triplet excita- 
tion spectrum in ID , leading to a Kondo spin liquid state for 
any finite values of coupling strength t/ J (with t as hopping 
and J as exchange) ; in 2D and 3D cubic lattices the spin gaps 
are found to vanish continuously around {t/ J)c ~ 0.70 and 
{t/J)c ~ 0.38, respectively, where quantum phase transitions 
occur and the Kondo spin liquid state changes into an antifer- 
romagnetically long-range ordered state. These results are in 
agreement with variational Monte Carlo, higher-order series 
expansion, and recent quantum Monte Carlo calculations for 
the symmetric Kondo lattice model. 

PACS numberes: 71.27.+a, 71.10.Hf, 75.10.Jm, 75.30. Hx 

Since the discovery of the class of stoichiometric insu- 
lating compounds known as Kondo insulators 0, there 
has been revived interest in the symmetric Kondo lattice 
Hamiltonian, 



(1) 



as a model of concentrated magnetic impurity spins cou- 
pled to conduction electrons. One of the important is- 
sues is the interplay between the Kondo screening and 
the magnetic interactions among localized spins mediated 
by the conduction electrons. The former effect favors a 
nonmagnetic Kondo spin liquid (singlet) state, while the 
latter interactions tend to stabilize an antiferromagnet- 
ically (AF) long-range ordered state. The character of 
such a transition between these two distinct phases has 
been a long standing issue since it was first pointed out by 
Doniach There have been a lot of investigations for 
the symmetric one-dimensional model, showing that its 
ground state is a disordered Kondo spin liquid state for 
any finite values of the coupling strength t/J [^. For two 
and three dimensional models, however, various approx- 
imate approaches such as variational Monte Carlo calcu- 
lations higher-order series expansions quantum 
Monte Carlo simulations and infinite dimensional 

calculations suggest that the Kondo spin liquid state 
may change into an AF long-range ordered state at cer- 
tain value of the coupling strength at low temperatures. 



Since there are a lot of difficulties in directly attack- 
ing the symmetric Kondo lattice model even in the ID 
case, a simplified version called Kondo necklace model 
was introduced by Doniach |2| , 
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where both and Si are spin 1/2 Pauli operators, denot- 
ing the conduction electron spin and impurity spin oper- 
ators, respectively, and means summation over the 
nearest neighbor conduction electron sites. Actually this 
simplified model is meaningful in general D dimensions 
{D = 1, 2, 3) in its own right. Due to the suppression of 
charge fluctuations in the symmetric model, the charge 
degrees of freedom are frozen out, so the first term of 
Eq.(2) represents the spin degrees of freedom imitating 
the propagation of the conduction electrons. This can be 
clearly seen in the ID case, where the first term is equiv- 
alent after a Jordan- Wigner transformation to a band of 
spinless fermions, which interact with localized spins via 
an AF spin-spin exchange coupling . 

Although the simplified model has only U(l) spin sym- 
metry, lower than SU(2) for the Kondo lattice model, 
the essential feature of these two models is kept. Thus, 
one would expect that the main physical properties of 
the original symmetric Kondo lattice model should be 
maintained in the Kondo necklace model. However, 
most of approaches used to treat the ID Kondo neck- 
lace model, including the variational mean-field calcula- 
tion approximate real-space renormalization group 
theory Q, and recent finite size scaling analysis [ p^ , 
have found a finite critical value of coupling strength 
{J/t)c ~ 0.24 — 0.38, below which an AF quasi-long-range 
order state appears, in contrast to Jc — 0, the result 
of quantum Monte Carlo simulation for the ID Kondo 
necklace model and the numerical result for the ID 
symmetric Kondo lattice model [|l^,|l^ . It is thus contro- 
versial whether the simplified spin model can be used to 
approximate the original symmetric Kondo lattice model. 
In this paper, we try to resolve this issue, starting from 
the Kondo necklace model, using the Kondo spin singlet 
and triplet representation, to reproduce correct ground 
states of the symmetric Kondo lattice model. In the ID 
case, the system is found to be in a Kondo spin liquid 
state with a finite spin gap for any finite t/J, while on 
2D and 3D cubic lattices a quantum phase transition 
occurs around {t/J)c ^ 0.70 and {t/J)c ^ 0.38, respec- 
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tively, where the Kondo spin hquid state changes into an 
AF long-range ordered state, in excellent agreement with 
the variational Monte Carlo calculation |^ , higher-order 
series expansion Q], and recent quantum Monte Carlo 
simulation on the corresponding symmetric Kondo 
lattice model. 

Our starting point is the strong coupling limit t = 0, 
where the lowest energy state of the model Hamiltonian 
Eq.(2) reduces to a sum over contributions from inde- 
pendent local Kondo spin singlet states at each lattice 
site. When t ^ Q, interactions between these indepen- 
dent local Kondo spin singlets are switched on. It will 
be seen later that this leads to very reasonable results 
even for t > J, which is of interest here. Usually, for two 
s = 1/2 spins Ti and Si placed on a lattice site, the local 
Hilbert space is spanned by four states consisting of one 
singlet and three triplet states defined as being created 
out of the vacuum |0) by the singlet and triplet creation 
operators: \s) = s^|0) and \ta) — tl^\0), {a = x, y, z). A 
representation of the impurity spins and conduction elec- 
tron spins in terms of these singlet and triplet operators 
is given by 

where a, [3, and 7 represent components along the x, ?/, 
and z axes, respectively, and e is the antisymmetric Levi- 
Civita tensor. This type of spin representation in terms of 
singlet and triplet (bond) operators was first proposed by 
Sachdev and Bhatt to study the properties of dimerized 
phases |l4j and then it was successfully used to consider 
the spin ladders |15| and s = 1 antiferromagnetic Heisen- 
berg spin chains pq |. As shown later, this representation 
faithfully describes the low temperature physics in the 
symmetric Kondo lattice model. In order to restrict the 
physical states to either singlets or triplets, a local con- 
straint is introduced: s\sn + a^n," = 1- Taking the 

singlet and triple operators at each site to satisfy bosonic 



commutation relations: 



and [s„,tj^ ^ = 0, the SU(2) algebra of the spins t„ and 
Sn can be reproduced. 
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[S'n,Q,T„,^j] = 0, S^=T^ = -. (4) 



Substituting the operator representation of the impurity 
and conduction electron spins, we obtain the following 
form of the model Hamiltonian, 

H = Hq + Hi + H2 + -ffa, 
J 



i a 

+ + - 1), 



+ U,yt]^y ] + h.C. 



^s(t 

^i.x 



t 

A- 



f ( f A- f - 



h.C. 



where a site-dependent chemical potential has been 
introduced to impose the local constraint. Here the local 
spin triplet states are split into two parallel spin states 
with TOs = ±1 and an anti-parallel spin state with rus = 
0. Hi describes the couplings between the singlet state 
and the parallel spin triplet states, while H2 corresponds 
to the couplings of the parallel spin and the anti-parallel 
spin triplet states. Hy, describes an interaction of one 
singlet boson and three different components of triplet 
bosons. 

The above Hamiltonian can be solved by a mean field 
decoupling of the quartic terms. It yields an effective 
Hamiltonian H^f with only quadratic operators. We 

take {s\) ~ {si) ~ s, which corresponds to a conden- 
sation of the local Kondo spin singlets on each site in 
accordance with the configuration of the ground state in 
the strong coupling limit, and the local chemical poten- 
tial is replaced by a global one. We will consider here 
only the terms Hq and Hi , as it can be shown that inclu- 
sion of H2 changes the results only slightly |15 1^ and 
all the decouplings of H3 identically vanish within the 
present mean field theory. After performing a Fourier 
transformation of the boson operators, ti^a — ^k.a 

k 

g-ik-Fi^ the mean field effective Hamiltonian is given by 
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mf 
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-J S^ + fiS^ 



J2 [Ak4yk,/3 + Ak (4 



k,/3— a;,y 

with Ak = + |is^A(k), Ak 



jts'^X(k), and 



M^) — cos ka. The lattice spacing has been taken to 

a=l 

be unity. This mean field Hamiltonian can be diago- 
nalized by a Bogoliubov transformation into new boson 
operators: ik,/3 = Wktk,/3 + "^^k^Lk/?; where the coeffi- 
cients Uk and Wk are even functions of k, and are de- 



termined to be: I 
2wkWk = sinh 29]^ 



v? = cosh26'k = 



2Ak 



VAj-(2Ak)2 



Then we obtain 



Hmf = t^oX!*L*k,z 



E 

\<.,f3—x,y 



Ea 



and 



(7) 
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where ujo = ( j + is the dispersionless energy level 
of the anti-parallel spin triplet excited state, Wk — 
— (2Ak)^ corresponds to the excitation spectrum 
of the parallel spin triplet excited states, and the ground 
state energy of the system is Eg = TV (— | J + /is^ — + 
^(cjk^Ak). By minimizing the ground state energy 

k 

with respect to /i and s, we derive the following saddle- 
point equations, 



N ^ 



A. 
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When a dimensionless parameter d 



is intro- 



duced, a self-consistent equation for d can be obtained 



can also be calculated from the saddle-point equations. 
The minimum spin gap appears in the parallel spin triplet 
excitation at k = (tt, tt): A^p — (i + j) %/! — Id, dis- 
played in Fig. 2. The most important feature here is that 
as the coupling parameter tj 3 increases, the drop of the 
spin gap in the small values oitjj continues down to 
the point {t/ J),. « 0.70 where the spin gap actually van- 
ishes. The critical coupling {t/J)c ~ 0.70 corresponds 
to a quantum critical point for a phase transition from 
the quantum disordered Kondo spin liquid to a magneti- 
cally long-range ordered state. Surprisingly, the location 
of the critical point for the 2D Kondo necklace model is 
precisely the value obtained from the variational Monte 
Carlo calculation the higher-order series expansion 
Q, and recent quantum Monte Carlo simulation for 
the 2D symmetric Kondo lattice model. When a simi- 
lar calculation is carried out in the 3D Kondo necklace 
model, the minimum spin gap appears at k = (tt, tt, tt) 
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and Asp = J {{ + j) Vl - 3d, shown in Fig.3. As t/J 



(9) 



to determine the variational parameters s and jj, and the 
spin triplet excitation spectra: wo = <^(|:+j) and ujk = 
J (I; -|- j) yjl + dA(k). There is a minimum spin gap in 
the parallel spin triplet spectrum at the AF reciprocal 
vector momentum k = Q: A^p = J (i -|- ^) ^1 — Zd/2, 
where Z is the total number of the nearest neighbors on 
the cubic lattice. 

In the ID case, we first numerically calculate the pa- 
rameters d, s^, and ji/J for a range of the coupling 
strength < t/J < 5, and the minimum spin gap 
Asp — J {j + ^) — d is evaluated in the range of 
< t/J < 5, which has been delineated in Fig.l. The 
dispersive band can also be parameterized by a spin den- 



grows, the spin gap decreases and exhibits a critical value 
{t/J)c « 0.38, where the spin gap disappears completely, 
showing a quantum phase transition from the quantum 
disordered Kondo spin liquid to a magnetically long- 
range ordered state as well. This transition point is in 
the same range as the higher-order series expansion Q for 
the 3D symmetric Kondo lattice model: {t/J)c ~ 0.50. 

Moreover, the present mean field theory can also be 
applied to the magnetically long-range ordered phase in 
the 2D and 3D Kondo necklace models. If we assume 
that, not only the local Kondo spin singlets (s bosons) 
condenses, one of the local Kondo spin triplets (tk.^ 
bosons) condenses as well on the AF reciprocal vector 
'•k X = V^^f^k.Q + '7k,K, corresponding to /jxmi? the orien- 
tation of the localized spins along x- direction, it will lead 
to another mean field effective Hamiltonian 



sity wave with a velocity given Vs = J + j') J ^. A H'„^j = E'^ + ujo^^ tLz^Kz + E ^k(^k,y^i^,v + '^Lx^K^) , 



linear drop of the spin gap is seen for small values of t/J. 
As t/J gets larger, the spin gap deviates considerably 
from the linear behavior and there is no indication at all 
suggesting a critical value for t/J where the gap would 
vanish. Since the excitation spectra are real and positive 
everywhere in the Brillouin zone, the system will be in a 
quantum disordered — Kondo spin liquid state for finite 
values of the coupling strength t/J, and the spin-spin cor- 
relation function decays exponentially at large distances 
with a correlation length ^ — 4^. This is indeed consis- 
tent with both quantum Monte Carlo simulations for the 
ID Kondo necklace [?] and numerical results for the ID 
symmetric Kondo lattice model |^,^. Thus, starting 
from the limit t/J = with localized Kondo spin singlets 
on each site, we see that any finite coupling strength delo- 
calizes the local Kondo singlets, reducing the magnitude 
of the gap but not closing it completely. 

Having secured the correct ground state for the ID 
symmetric Kondo lattice model, we now turn to two and 
three dimensional "Kondo necklace" models on a cubic 
lattice. In 2D, the variational parameters d, s^, and /i/ J 



E'g^N 



-Js + /is — /i 



4 " V 4 

(wk - Ak) , 



J , 1 „ ,_2 \ t2 

/i Zts I t 



(10) 



where Wk has the same form as in the Kondo spin liquid 
phase, and rj^^ ^ and ^k,x are the transverse spin triplet 
excitation mode. When the order parameter t is nonzero, 
the saddle point equation for t yields /i = \Ztl^ — j, 
which makes the parallel spin triplet excitation spectrum 
gapless: Wk — ^ZfE^^Jl + 2A(k)/Z. The ground state 
corresponds to a magnetically long-range ordering state 
with a maximum momentum q = Q, and the mean field 
t represents the AF order parameter. It has been sug- 
gested that a very appealing physical picture of forming 
AF long-range order in the Kondo necklace or the sym- 
metric Kondo lattice models: when t/J is small, the con- 
duction electron spins are locked and the impurity spins 
are screened completely, and the ground state is a prod- 
uct of the local Kondo spin singlets - quantum disordered 
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phase g. Ast/J becomes larger and larger, the conduc- 
tion electrons (the spin degrees of freedom) have more 
possibility to propagate to the nearest neighbor sites, and 
the localized magnetic impurity spins is only partially 
screened (s ^ 0), then the remaining part of the mag- 
netic impurities on different lattice sites start to develop 
long-range correlations {t ^ 0) mediated by the conduc- 
tion electron spins |p^ . Such a magnetically long-range 
ordered state might be related to the ground states of the 
U-based heavy fermion compounds (URu2Si2 and UPta) 
with a very small magnitude of induced staggered mag- 
netic moments. In order to determine the parameters s 
and t, we minimize the ground state energy, derive the 
saddle point equations, and finally obtain 



J 

J 



^5]v/l + 2A(k)/Z, 



2N 



,/TT2: 



y/l + 2X{k)/Z 



The AF order parameter is defined by nis 
to the following expressions: 



rus = ^/ (0.35712 - ^) (0.52095 -I- ^) 
= \/ (0.44234 - ^)(0.51263+ ^) 



(11) 
s t, leading 

for 2D; 
for 3D. 



These results have also been displayed in Fig. 2 and Fig. 3, 
respectively. In Fig. 2, our results are also compared with 
the numerical results for the spin gap and staggered mo- 
ment of the magnetic impurity spins in the recent quan- 
tum Monte Carlo simulation on the 2D symmetric Kondo 
lattice model at zero temperature |^ . 

In summary, we have presented a mean field theory 
for the Kondo necklace model in ID, 2D and 3D and 
have obtained their correct ground states corresponding 
to the respective Kondo lattice model. A long standing 
controversial issue has been thus resolved regarding the 
relationship between these two models. As far as the 
spin part of the ground state properties is concerned, the 
Kondo necklace model can reproduce the correct phase 
diagrams of the symmetric Kondo lattice model at zero 
temperature. 
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Figure Captions 

Fig.l. The variation of the spin gap upon increasing 
of the coupling parameter t/ J of the ID model at T = 0. 

Fig. 2. The spin gap and the staggered magnetic mo- 
ment at zero temperature of the 2D Kondo necklace 
model (bold line) in comparison with results of recent 
quantum Monte Carlo simulation ^ for the 2D Kondo 
lattice model. 

Fig. 3. The spin gap and the staggered magnetic mo- 
ment at zero temperature for 3D Kondo necklace model. 
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